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Abstract 

Two models of binary fragmentation are introduced in which a time depen- 
dent transition size produces two regions of fragment sizes above and below 
the transition size. In the models we consider a fixed rate of fragmentation 
for the largest fragment and two different rates of fragmentation for the two 
regions of sizes above and below the transition size. The models are solved 
exactly in the long time limit to reveal stable time-invariant solutions for 
the fragment size distributions. A rate of fragmentation proportional to the 
inverse of fragment size in the smaller size region produces a power law dis- 
tribution in that region. A rate of fragmentation combined of two terms, one 
proportional to the inverse of the fragment size and the other proportional to 
a logarithmic function of the fragment size, in the larger size region produces 
a log-normal distribution in that region. Special cases of the models with no 
fragmentation for the smaller fragments are also considered. The similarities 
between the stable distributions in our models and power law log-normal dis- 
tributions from experimental work on shock fragmentation of long thin glass 
rods and rupture of mercury droplets are investigated. 
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1 Introduction 



Fragmentation is an irreversible kinetic phenomenon which occurs in many 
physical and chemical processes. Because of broad range of applications, 
many recent studies have been carried out to investigate the kinetics of the 
processes by introducing simple fragmentation models [HHKj. In [|l|-|3| Ziff 
and MaGrady presented a model of fragmentation in which the rate of break 
up depends on the size of the fragments. A general discussion of the kinetics 
of continuous fragmentation processes was given by Cheng and Redner in || . 
In [7],H Krapivsky and Ben-Nairn studied the kinetics of random fragmenta- 
tion of multidimensional objects. In || three models of binary fragmentation 
were investigated analytically in which at each time step, the largest frag- 
ment in the system is broken with some externally tuneable probability. In 
the second model of that paper, at each time step, there is a fixed rate of 
fragmentation for the largest fragment and a rate proportional to the inverse 
of the fragment size for all other fragments smaller than the largest one. 
The model was solved exactly in the long time limit to reveal a power law 
distribution with an exponent, which depends on the precise details of the 
fragmentation process. 

Surprisingly, there is a similarity between the result of the stable distri- 
bution from the second model of ||, with some experimental and analytical 
results of fragment size and mass distributions with power law forms in the 
shock fragmentation |11|-20||. In JT(| a generalization of the second model 



of || was carried out. Some models of binary fragmentation were introduced 
which revealed composite power law distributions for the mass and size of 
the fragments flOfl . 

Another behaviour of distribution which is observed in many different 
processes has a log-normal form |T^p^|. In a comprehensive study in 
the log-normal distribution is given. The application of these distributions in 
many different areas of science is discussed in , such as economics, biology, 
astronomy, philology, small particle statistics and physical and industrial 
processes. The log-normal behaviour is known to be able to describe the size 
distribution well in a wide variety of geological situations [p4j] , such as that of 
rocks in a boulder field. In [18] an analytical and numerical study was carried 
out on the fragmentation of long thin glass rods. Some stochastic models for 
one-dimensional brittle fracture were introduced. The models successively 
extended to describe cascade fracture of long thin glass rods. They are 
found to give the fragment size distribution close to a log-normal one. In 
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terms of rupture in liquids, experiments in rupture of mercury drops were 
performed ||19|| . The result of size distribution for the cumulative number 



of droplets showed a log-normal distribution at the small falling heights. 
A transition of the distribution from log-normal to scaling behaviour was 
observed as the falling height was increased. In |2(J an experimental work 
on the shock fragmentation of long thin glass rods was carried out by Ishii 
and Matsushita. The results of fragment size and mass distributions at small 
falling heights showed a log-normal form for larger fragments and a power 
law form for smaller fragments. The crossover was seen to be at length scales 
around the rod diameter. The mass and size distributions for the larger 
fragments showed a power law form as the falling height was increased. 

This paper is organised as follows. In sections ^| and ^| we introduce two 
new models, which are generalisations of the second model of |J and the mod- 
els of fragmentation in |TIJ. In the models we introduce a time dependent 
transition size which produces two regions of fragment sizes with different 
rates of fragmentation above and below the transition size. A rate of frag- 
mentation proportional to the inverse of the fragment size in the smaller size 
region produces a power law distribution in that region. A rate of fragmenta- 
tion combined of two terms, one proportional to the inverse of the fragment 
size and the other proportional to a logarithmic function of the fragment size 
yield a log-normal distribution in the larger size region. Finally in section [| 
we summarise our conclusions and discuss the application of our results to 
experimental work on the shock fragmentation of long thin glass rods GO 
and the rupture of mercury drops ||19|| . 

For the rest of this introduction we give a very brief summary of previous 
work on models of binary fragmentation. In these models, the density of 
fragments of size y at time t, n(y,t), evolves according to 



n(y,t + 5t) = n(y,t) — 5tn(y,t) / R(z,y—z)dz + 2 5t / R(y, z—y) n(z,t)dz, 

Jo Jy 

(1) 

where R(y, z) is the intrinsic rate that a fragment of size (y + z) breaks 
into fragments of size y and z. The first term on the right hand side is 
the contribution from those fragments not chosen for the fragmentation in 
time (t, t + 5t). The second term represents the decrease in the number 
of fragments of size y by the fragmentation into fragments of size z and 
(y — z)(y > z). The third term represents the increase in the number of 
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particles of size y due to the fragmentation of particles of size z(> y), such 
that one of the products is of size y. 

The upper limit of the second integral on the right-hand side is usually 
set to some fixed value greater than the size of the largest particle in the 
system. This value is usually 1 or oo depending on the details of the model. 
To connect with the analysis of our models in the following sections, we have 
set the upper limit to L(t), the size of the largest particle at time t. As 
n(y,t) = for all y > L(t), this choice has no effect on our results. 

We now turn to a detailed investigation of our models. 



2 Model A 

Firstly, we introduce a time dependent transition size, y m (t), between zero 
and the size of the largest particle in the system at time t, L(t). The transi- 
tion size produces two regions of fragment sizes with different rates of frag- 
mentation above and below the transition size. At each time step there is 
a fixed probability, p±, for the fragmentation of the largest particle in the 
system. A particle of size smaller than the transition size breaks with a 
probability p 2 and a rate proportional to the inverse of the fragment size, 
at each time step. Any other particle of size larger than the transition size 
can be fragmented, at each time step, either with a probability p 3 and a rate 
proportional to the inverse of the fragment size, or with a probability p 4 and 

a rate proportional to — — In [ ^-A- ) • Here y(t) is the size of the fragment at 

y(t) \y(t)J 

time t and y(t) is a size for which the logarithmic rate is zero at time t. In the 
model we have P1+P2 + P3+P4 = 1- We also assume that y m (t) > y(t), such 
that for all fragments of sizes larger than the transition size, the logarithmic 
term in the rate of fragmentation is positive. In order to make the model 
tractable, we will later choose to have y m (t) proportional to L(t). 

In the model the distribution function of daughter fragments is uniform. 
Therefore the rate of breaking of a particle into two smaller pieces is indepen- 
dent of the daughter fragment sizes. With probability p 2 every particle of size 
smaller than y m (t) is equally likely to be chosen for the fragmentation. With 
probabilities p 3 and p 4 , every particle of size larger than the transition size is 
equally likely to be selected for the fragmentation with the correspondence 
rate. We assume that in time (t, t + St) the largest size changes from L(t) to 
L(t) — p\SL and the transition size changes from y m (t) to y m (t) — Sy m . 
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Then, the fragmentation rate for this model is given by 



5tR(y,z) 



I'l T77V Kv + z ~ L (t)) + P2 0(y m {f) - (y + z)) 



+ 



Lit) 

St St , 

P3 — ; h Pi — ; — In 

y+z y+z 



y + z 

' y + z 



6{{y + z)-y m {t)). 



(2) 



The 5-function in the first term on the right-hand side ensures that only the 

SL 

largest particle is fragmented and -— - is the probability of placing the cut 

L[t) 

in the largest particle in a particular infinitesimal length SL. The second 
term on the right-hand side of @ represents the rate of fragmentation for 
the particles of sizes smaller than y m (t) in time (t,t + St). The remaining 
terms on the right hand side of (0) represent the rate of fragmentation for 
the fragments of sizes between y m (t) and L(t) in time (t, t + St). Inserting the 
rate (0) into the kinetic equation (0) for the fragments in both size regions 
smaller and larger than y m (t), yields 



m(y,t + St) = ni{y,t) - p 2 Stn 1 (y,t) + 2p x 



SL 
W) 



n 2 (L(t),t) 



c f ym(t> , s dz „ f L W . . dz 
+ 2p 2 St n x (z,t) — +2p 3 St / n 2 (z,t) — 

rL(t) ( Z \ dz 

+ 2p 4i St / n 2 (z, t) In 
Jy m (t) 



(3) 



for y < y m (t) and 



n 2 (y,t + St) = n 2 (y,t) - p 3 Stn 2 (y,t) -p A St\n 



n 2 (y,t) 



SL r L (t) dz 

+ 2 Pl —--n 2 {L{t),t) + 2p 3 St ^ n 2 {z,t) — 



f L{t) 

+ 2p 4 St / n 2 (z, t) In 
Jy 



z \ dz 



Mt)J z 



(4) 



for y > y m (t). Here rii(y, t) and n 2 (y, t) are the densities of fragments of size 
y at time t which is smaller and larger than y m {t), respectively. The second 
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term on the right hand side of (Q) and the second and third terms on the 
right hand side of (f|) represent the decrease in the number of particles of 
size y from fragmentation into smaller particles, with probabilities p 2 , P3 and 
£>4, respectively. The next term on the right hand sides of ([|,fD describes the 
gain in the number of particles of size y from the fragmentation of the largest 
fragment of size L(t). The remaining terms on the right hand sides of (|3]|D 
are the increase in the number of particles of size y from the fragmentation 
of all particles larger than y and smaller than L(t). 

Now we define densities of fragments of length y at time t, which are 
normalised to a positive constant value A for two regions of sizes smaller and 
larger than y m {t), 

^vA = A 1 ^ ) "-^jm («) 

/ rii(z,t)dz + / n 2 (z,t)dz 
Jo Jym(t) 

for y < y m (t) and 

( + \ a n 2 (y,t) 

^>*)= A 7« JW) ( 6 ) 

/ ni(z,t)dz + / n 2 (z,t)dz 
Jo Jym(t) 

for y > y m {t). In time t — >• t + 5t as L(t) — > L(t) — p\5L and y m {t) — > 
(t) -5y m they yield 



gi{y,t + 5t) = A 



/ rii(z,t + 5t)dz + / n 2 (z } t + 5t)dz 

JO Jy m {t)-Sy m (t) 

(7) 

for y < y m (t) and 



/ . i tja a n 2 {y,t + St) 

g 2 {y,t + dt) = A 



rym(t)-dy m (t) rL[t)-p 1 dL 

/ ni(z, t + 5t) dz + / n 2 (z,t + 5t)dz 

JO Jy m {t)-Sy m (t) 

(8) 

for y > y m (t). Inserting ([3]), (f|) in (^), (|), respectively, gives us 
g 1 {y,t + 8t) = 



SL rym(t) 
2 Pl —^-g 2 (L(t),t) + 2p 2 5t I 9l (z,t) 



[L(t) d z r-i^yt) 

+ 2p 3 St g 2 (z,t) h 2p A St / g 2 (z,t)^, , 

Jy m (t) z J ym (t) \y(t) J z 

l- P2 St- ^g 2 (L(t),t)5L -^Stj o gi(z,t)di 



dz 



dz 

y z 

L(t) 



In 



z \ dz 



+ gi(y,t) 



A Jy m (t) 



fL(t) p A rL{t) ( Z \ 

/ g 2 (z,t)dz -St g 2 (z, t) In — — dz 

Jym(t) A Jy m (t) \y[t)J 



(9) 



for y < y m (t) and 



g 2 (y,t + St) = 

SL /- L W dz 

2 Pl —^g 2 (L(t),t) + 2p 3 8t ^ g 2 {z,t) — 

rL(t) ( Z \ dz 

+ 2p±StJ^ g 2 (z,t)\n ( — ] — + g 2 {y,t) 



— Pi St In 



y 



- ^g 2 (L(t),t)SL-^St 
y(t) J A A Jo 



1 - ps St 

Vm(t) 



gi(z,t)dz 



P3 f L ^ P4 f L ^ 

— St g 2 (z,t)dz - —St g 2 (z,t)\n 

A Jy m (t) A Jy m (t) 



z 



Mt). 



\dz 



(10) 



for y > y m (t). In order to solve (|) and (pTOj ) we introduce two functions 
Fi(x, t) and F 2 (x,t) defined by 



F 1 (x,t) = L{t)g 1 (xL{t),t) 



(11) 



for x < x, m and 



F 2 {x,t) = L{t)g 2 {xL{t),t) 



(12) 



for x > x m . Here we have defined dimensionless and time independent 



variables x 



y(t) 



X ■■>■ 



y m (t) 



and x 



where x changes in the range 



L(tY m L(t) L(t) : 
[0, 1]. x m and x are fixed values and restricted to the same range. Using (fi- 
ll!?!) we can obtain two partial differential equations for F 1 (x,t) and F 2 (x,t) 
as 
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dF^t) 
dt 



— vx 



dF 1 (x,t) 
dx 



+ 2uF 2 (l,t) + 2p 2 I F 1 (z,t)— +2p 3 I F 2 (x,t) 



dx 



x 



+ 2p 4 I 1 F 2 (x,t)\n(%\— - F 1 (x,t) 

Jx m \X / X 



v + P2 + jF 2 (l,t) 



+ 



P2 

A 



F 1 (x,t)dx + 



A 



F 2 (x,t)dx + ^ 

Xm A 



F 2 (x,t) In ( - ) dx 



x 



(13) 



for x < x m and 



dF 2 (x,t) _ 



dt 



V + P3 



at 

dF 2 (x,t) ^ . . r 1 ^ , , dz 

-vx 2 . V ' ' + 2i/F 2 (l,t) + 2p 3 / F 2 {z,t) — 

OX J x Z 

+ 2 Pi I' F 2 (z,t)\n(^)-- F 2 (x,t) 

Jx \X) Z 

{ X\ V Vo f Xm 

+ p 4 In + - F 2 (l, *) + ^ / Fi(x, t)dx 
^ f 1 F 2 (x,t)dx + ^ ^ F 2 (x,t)\n 

to be a fixed value, and assumed 



+ 



dx 



(14) 



SL 



for x > x m . Here we have set Pi-r— - 

L(t)8t 

that pi > 0. i/ is equal to the rate of fragmentation for the largest fragment, 
and determines the relationship between the real time and the length of the 
largest fragment at time t, L(t). In the limit t — > oo, we assume that Fi(x, t) 
and F 2 (x,t) evolve to time-independent quantities, 

F 1 (x) = \im t _ >QO F 1 (x,t) (15) 

for x < x m and 



F 2 (x) = lim^oo F 2 (x,t) 



(16) 
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f q j-i. j- . dF^t) dF 2 (x,t) .... 
tor x > x m . bo that, as t — > oo then ^- > and — > 0. lo 



dt 



dt 



solve (|13|Jl^) in the long time limit using (|T5l , [T6|) , we anticipate solutions as 

F 1 (x)=Bx- a (17) 



for x < Xm and 



F 2 (x) =C 



exp | - [in (§) 



x 



(18) 



for x > x m . Substituting (0,0) in (|I3|]I4[) in the long time limit, using 
( |T3|JTE| ) give us 



P2 

a = — , 

P3 







Pi 
2p 3 



(19) 
(20) 



and 



(21) 



B and C in fllTD , (|lq) , respectively, are constant values and can be obtained 
using (|l"3l , |l4l) in the long time limit and the continuity condition at x m , 



FAx r . 



Fo(x r . 



(22) 



These give 




9 



where the second and third terms in the denominator of ( |2"3] ) are a proportion 
of the normal probability distribution function defined as 



P(x) 

We have also assumed that 



t2 ' 2 dt. 



(25) 



< a = — <1. (26) 
Ps 

Solutions flTFyiBD using (p^, p0| , p^ -p6|) satisfy the normalisation relation 

/ F 1 (x)dx + / F 2 (x)dx = A. (27) 

J0 J x m 

If A = 1, the normalisation is equal to 1. For a choice of A as 



A 



Pi 



P3 



■ cxp 



hi 



X 



2vrp 3 p 3 - P2 
the solutions usin S HHHH yield 



P3 



P I J— In I — 



P3 



(28) 



*i(z) 



27rp 3 



exp < 




In 


•Em \ 
\ X ) _ 






) 



— E2. 

X p 3 



(29) 



for x < Xm and 



Fo(x) 



' Pa 
2np 3 



exp 



Pi 
2ps 



In 



1 2 



(30) 



for a; > x m . We see that the stable distribution has a power law behaviour 
in the smaller region, x < x m . For a log- normal distribution, the probability 
of finding a fragment of size between x and x + dx is given by n(x)dx, where 



n(x) 



cxp 



{-[In {x/x)} 2 /2a 2 } 



(27ia 



2^/2 



(31) 
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where x is the mean and a is the dispersion of the distribution. Comparing 
p0| ) with ( |3~1~D shows that the stable distribution F 2 (x) in the larger region 
exhibits a log-normal behaviour with a mean and a dispersion of the distri- 
bution equal to x and a = J—, respectively. To have a positive value for A 

^ Pa 

in ( p8|) and therefore physically meaningful distributions, we require P3 > £>2- 
This requirement is satisfied in (Pfi|). 

Now we consider some special cases of the model. 

1. There is no fragmentation for the fragments smaller than y m (t), i.e. 
P2 = and pi + P3 + pi = 1. The stable distribution then from fl2~9[ ) 
yields 



/ Pi 
2np 3 



exp 



P4 
2^3 



In 



X 



x 



X : 



(32) 



for x < x m and is the same as ( p0|) for x > x m . We see that the 
distribution is a constant value in the smaller region and has a log- 
normal form in the larger region. 



2. As y m (t) = y(t) and therefore x m = x, the distributions (^,^3) give 



Fi(x) 



1 Pa 
2ttj9 3 



Xr, 



P3 



X p 3 



(33) 



for x < Xm and 



Fo(x) 



I Pa 
2vrp3 



exp< - 



Pa 
2p 3 



In 



x 



X i 



X 



(34) 



for x > x m . Again the distribution reveals a power law form in the 
smaller region and has a log-normal behaviour in the larger region. 

P3 



3. As p4 — > and therefore a 2 



oo, the distribution (BUD in the 



larger region yield F 2 (x) 



Pa 



~ x 



On the other hand as p^ — > the 



rate of fragmentation (|]) for the fragments of sizes larger than y TO (i) 
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is only proportional to the inverse of the fragment size. Consequently, 
at each time step either the largest fragment breaks with probability 
Pi, or a fragment smaller than y m (t) is chosen for fragmentation with 
probability p 2 , or a fragment larger than y m (t) is selected with proba- 
bility P3, where pi + p 2 + P3 = 1- The stable distributions then exhibit 
a composite power law behaviour with power laws of exponents equal 

to — — and — — in the smaller and larger region, respectively In our 
v v 

case since from (Ell) we have v = p 3 , the exponents reduce to and 

—1 in the two regions. The power law distribution in the smaller region 
is similar to that from (|l7],|rj|) and in the larger region has an exponent 
of —1 as expected. The latter is consistent with one of interesting prop- 
erties of the log-normal distribution where as the variance a 2 is pretty 
large, then the distribution can be approximated by a power law form 
with an exponent equal to —1, pX 



3 Model B 



In this model again a time dependent transition size, y m (t), produces two 
regions of fragment sizes with different rates of fragmentation above and 
below the transition size. At each time step either the largest fragment in 
the system breaks with a fixed probability, pi, or a fragment smaller than 
the transition size is chosen for the fragmentation with probability P2 and a 
rate proportional to the inverse of fragment size, or another fragment larger 
than y m (t) is selected with probability (1 — pi — p 2 ) and a rate proportional 



to 



St 

W) 



1 + A In 



positive and constant va" 



Here y(t) is the size of fragment at time t, A is a 



ue and y(t) is a size for which the logarithmic term 
is zero at time t. In this model we assume that y m (t) < y(t). For a fragment 
of size y > y(t) the logarithmic term in the rate of fragmentation is always 
positive, but a fragment of size y m (t) < V < y(t) gives a negative value for 
the logarithmic term. In order to have a positive rate of fragmentation for 



all particles larger than y m (t), we require 
y(t) > y(t) e" 1//A . If we choose y m (t) to be 



1 + Aln 



Mt). 



> which gives 



y m (t) = y{t) e 



l/A 



(35) 



12 



then the rate of fragmentation for all fragments of sizes larger than y m (t) is 
positive. 

Consequently, the rate of fragmentation for this model is given by 



6tR(y,z) 



Pi t^t % + * - Lit)) + p 2 6(y m (t) -(y + z)) 

L(t) y + z 



+ (I-P1-P2) 



St 



y + z 



1 + A In 



y + z 



6{{y + z)-y m {t)), 
(36) 



where y m (t) is defined as (|35|) . The first two terms on the right hand side 
of (|36|) are the same as those in (0) and have the same definitions as before. 
The last term represents the rate of fragmentation for the fragments of sizes 
larger than y m {t) in time (t,t + St). Using (|ll,|36l) for the both size regions 
give the kinetic equations 



n 1 (y,t + St) 



SL 



ni(y,t) - p 2 Stn 1 (y,t) + 2px—-rn 2 (L(t),t) 



dz 
z 



rvmw dz f^w 

+ 2p 2 St m{z,t)— +2(1- Pl -p 2 )5t 712(2, t) 
Jy z Jy m {t) 

r L (t) ( z \ dz 

+ 2\(l- Pl -p 2 )St n 2 (z,t)\n — - — (37) 
Jym(t) \y{t)J z 

for y < y m (t) and 



n 2 (y,t + St) = 

n 2 {y, t)-(l-px- p 2 ) St n 2 {y, t) - A (1 - p x - p 2 ) St In [jj^j n i{v, l ) 



SL f^ 1 ) dz 

+ 2p 1 —n 2 {L{t),t) + 2(l-pi -p 2 )8t / n 2 {z,t) — 

L(t) Jy Z 

r L (t) ( z \ dz 

+ 2 A (1 - Pl - p 2 ) St ^ n 2 (z, t) In — (38) 

for y > y m {t). Following the same procedure as that used in the model A, 
we obtain 
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dFi(x,t) 
dt 



— V x ■ 



9F 1 {x,t) 
dx 



+ 2uF 2 (l,t) + 2p 2 I F 1 (z,t) — 



+ 2 (1 - Pl - pa) F 2 (x, £) — + 2 A (1 - Pl - p 2 ) f F 2 (x, t) In (?) — 

-FiOM) 



f + P2 + -rF 2 (l,t) + ^ / F^x^dx 
A A Jo 



+ 



(1 -pi -p 2 ) 



A(l - pi - pg) 
4 



F 2 (x,t) In ( - ) dx 



(39) 



for x < Xm and 



dF 2 (x,t) 
dt 



V X 



dF ^ t>j +2vF 2 (l,t)+2(l- Pl - p 2 ) £ F 2 (z, *) 



rl /z\ dz 

+ 2A(l-pi-p 2 ) / F 2 («,t)ln F 2 0,t) 



Z/+ (1 - pi -p 2 ) 



+ A(l- Pl -p 2 ) In - + T F 2 (l,t) 



x 



X 



V 



A 



P2 

A 



Fi(x, t)dx 



+ 



(1 - pi - pa) 



F 2 (x, t)da; + 



A(l - pi -pg) 



F 2 (ar,t) ln( - ) dx 



x 



(40) 



for x > x m . Here x, x m , x and v have the same definitions as before and 
from Q55D we get 



x m = xe- 1 / A . (41) 

In the long time limit as t — > oo, we assume that Fi(x, t) and -F 2 (x, t) evolve 
to the same time-independent variables as (0) and (0), respectively. In 
order to solve (p9||40D in this limit we anticipate solutions as 



F 1 (x) =Gx- 



(42) 
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for x < x m and 



Fo(x) = H 



exp < — rj 



ln (l) 



x 



(43) 



for x > x m . Substituting ( f42l , [43"D in (|39| , |40|) in the long time limit, gives us 

P2 



7 



n 



l - pi - p 2 1 

A 
2 



and 



v 



l-pi- p 2 - 



(44) 
(45) 

(46) 



G and H in (0,^) are constant values and can be obtained using (|39],f|0D in 
the long time limit and the continuity condition at x m , (1221) . These give 



H 



and 



.4 



1 - Pi - P2 
l-pi- 2p 2 



1 

2A 



ex P - + V -r 




p(VAln(I 



P 



G = H 



eXPl "2A 



x e 



where we have assumed that 



< 7 



1/A\ 1— Pl— P2 



< 1. 



(47) 



(48) 



P2 



(49) 



1 - pi - p 2 

The second and third terms in the denominator of ( |4TD are a proportion 
of the normal probability distribution function with the same definition as 
Q25]). Solutions (^,^) using ( M - K5p satisfy the normalisation relation (P7|). 
If A — 1, the normalisation is equal to 1. For a choice of A as 
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A = \ — 



A I-P1-P2 



2vr 1 - Pl 
the solutions (H2 



2p 2 
using 



exp 



x 



P 



7x • (50) 



MfmHSI) yield 



for x < x e~ 1,/A and 




exp 



1 

2A 



l-p 1 -2p 2 
1/A\ 1-P1-P2 



X 



(51) 



x e 



^2(X) = W- 




exp 



A 
2 



lu ( - 



(52) 



for x > xe -1 / A . We see that the stable distribution has a power law be- 
haviour in the smaller region. Comparing (|52| ) with (31) exhibits that the 
stable distribution F 2 (x) in the larger region has a log-normal form with a 
mean of x and a dispersion of the distribution equal to 



a 



1 



(53) 



It is shown that the log-normal function ( ^2|) has a maximum at x = x e~ 
which we have assumed to be equal to x m . This shows that in this model 
the dimensionless transition size x m is chosen to be equal to the value for 
which the log-normal distribution in the larger region is maximum. To have a 
positive value for A in (|50|) and therefore physically meaningful distributions, 
it is required that (1 — p± — p 2 ) > P2- This inequality is satisfied from ([4"Up. 
Now we consider some special cases of the model B. 

1. When there is no fragmentation for the fragments smaller than y m {t), 
i.e. pi = 0, then the stable distribution from (|5T|) yields 



F 1 {x) = J±\e^ (54) 
V 2tt x 

for x < x m and is the same as (|52D for x > x m . The distribution is a 
constant value in the smaller region and has a log-normal behaviour in 
the larger region. 
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2. As A — > 0, and therefore from ( |5*3"D cr 2 — > oo, the situation is then simi- 
lar to the case 3 of model A and the log-normal distribution ([52] ) in the 
larger region can be approximated by a power law form with an expo- 
nent equal to —1, F2(x) ~ x~ 1 . On the other hand at this limit the rate 
of fragmentation ([36]) gives two different probabilities in the two differ- 
ent size regions with a rate proportional to the inverse of the fragment 
size in each region. The stable distribution then reveal a composite 

power law with power laws of exponents of — — and — — — in 

v v 

the smaller and larger size regions, respectively. Since in our case from 
(Hq) we have v = (1 —p% — P2), the exponents yield and 

1 - Pi ~ Pi 

— 1, respectively. The result in the smaller size region is equivalent to 



that from (42,53). In the larger size region, as expected, we get a power 
law with an exponent of —1 which is consistent with the property of a 
log-normal distribution with a very large variance. 

We see that all the results for the model B can be recovered from the 
model A by substituting p 3 by (1 — p\ — pz), p^ by A (1 — p\ — p<i) and using 
equation (|35|). 



4 Conclusions 

In this paper we have studied two models of binary fragmentation. In the 
models we introduced a time dependent transition size, y m {t), which pro- 
duced two regions of fragment sizes with different rates of fragmentation 
above and below the transition size. The rate of fragmentation at each time 
step is fixed for the largest fragments and is proportional to the inverse of 
fragment size for the fragments smaller than the transition size. We con- 
sidered a rate of fragmentation combined of two terms, one proportional to 
the inverse of the fragment size and the other proportional to a logarithmic 
function of the fragment size, for the fragments larger than the transition 
size at each time step. For a size equal to y(t) at time t, the logarithmic 
function is zero. In model A we assumed that the transition size, y m (t) is 
any arbitrary value for which y m (t) > y(t), whereas in the model B we chose 
a value as fl3"5] ) for the transition size. The models were then solved exactly in 
the long time limit to reveal stable time-invariant solutions for the distribu- 
tions. The results of these distributions for both models A and B exhibited 
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a power law form with an exponent between —1 and zero, in the smaller 
region. For a special choice of normalisation, the distributions exhibited log- 
normal behaviours in the larger regions. In this region they revealed one of 
interesting properties of a log-normal distribution for which as the variance 
of the log-normal is pretty large, then the distribution is approximated by a 
power law form with an exponent equal to —1. Special cases of the models 
with no fragmentation for the smaller fragments were also examined. The 
stable distributions then exhibited a constant value and a log-normal form 
in the smaller and larger regions, respectively. 

Now we investigate some experimental work on the fragmentation involv- 
ing log-normal distributions and the transition of the distributions from a 



log-normal form to a power law one and vice versa. In [EDI the experimental 
results of fragment size and mass distributions for the long thin glass rods 
with fixed lengths of 1500mm and diameters of 2mm were reported. The 
rods were dropped horizontally onto a flat hard floor from different heights. 
At lower falling heights (about lm drop), the distributions exhibited a log- 
normal form for larger fragments and a power law form for smaller ones. 
The crossover was seen to be at length scales around the rod diameter. This 
is due to the fact that fragments smaller than the rod diameter, undergo 
a three-dimensional fracture, whereas the larger fragments are produced by 
one-dimensional breaking. For a falling height of 1.20m the cumulative num- 
ber of fragments on their size were plotted. The data was fitted to a curve 
calculated by assuming that the size distribution could be described by a 
log-normal distribution. The fitting gave values of I = 30mm and a = 0.55 
for the mean and the dispersion of the log-normal distribution. The fitting 
is excellent for the fragments larger than l c = 7mm which is a length scale 
around the rod diameter. The size distribution of fragments smaller than l c 
seemed to have a power law form. As the falling height was increased, the 
cumulative number variation for larger fragments started to show a power 
law dependence on their size and mass. In [^] experiments in rupture of 
mercury droplets were performed. In the experiment at a height of h, the 
mercury droplets of about 2.0mm radius were fallen directly on a glass Petri 
dish. The results of the cumulative number of the drops versus the drop 
diameter were plotted. For small falling heights a log-normal behaviour was 
exhibited. The distribution showed a clear transition from a log-normal form 
to a scaled one as the falling height h was increased. 

Our models in the long time limit can be applied to the experimental 
systems when the fragmentation process is over and the system reaches a 
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stable and time-independent phase. In the models we assumed that the 
fragments have only one dimension, whereas in real physical systems the 
fragments have three dimensions. We can therefore only apply our models to 
the systems for which the fractures occur in one dimension. This assumption 
seems physically reasonable for a long thin glass rod |20| for which most of the 
fractures occur in the length of the rod and the cross-section remains almost 
fixed. We can also consider this assumption for the mercury droplets JT9J for 
which at each time step only the radius of the droplet changes and the shape 
of the fragments remains always spherical. 

Now we discuss the similarities between the results from our models and 
those from the experiments. The experimental results for the size distribution 
of long thin glass rods at low falling heights exhibited a power law and log- 
normal behaviour in the smaller and larger size regions, respectively |20 



We got the same behaviour for the distributions in both our models. The 
data from the experiments for the falling height of 1.20m revealed that l c < I, 
where l c is the transition size and I is the mean of the log-normal distribution. 
In the model B of our work also we have x m < x. We therefore expect model 
B of our work to describe the experimental results. In model B, the dispersion 
of the log-normal distribution can be obtained from ([5^) and the mean and 
the transition values, x and x m , respectively, are any arbitrary values which 
satisfy the relation ([H]). Then, by appropriate choices of x and A in the model 
B of our work we can match the mean and the dispersion of the log-normal 
distribution in the experimental results. Because of (ET|,F1) we should expect 



a relation as l c = I e~ a . The values of l c and I from experimental result do 
not satisfy this relation. As a result, although model B of our work does not 
match the exact relation between l c and I from experimental result, but it 
can give an explanation for the transition of distribution from a power law 
to a log-normal form on the shock fragmentation of long thin glass rods ||20|| . 
Also in the smaller region of the model B, the appropriate choices of p\ 

P2 

and P2 can match the exponent of power law equal to — from 

(l-p 1 -p 2 ) 

(|42]j4^) , with the corresponding exponent of power law distribution from the 
experimental results. This exponent needs to satisfy the condition (^). For 
high falling heights a transition of the distribution from log-normal to scaling 
law was observed in the experiments for the larger rod fragments [|20| and for 
the mercury droplets [0. At high heights therefore the distributions from 
the experiments exhibit a composite and a single power law form for the 
rod fragments and mercury droplets, respectively. The models in this paper 
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suggest a power law and a log-normal behaviour for the distribution of the 
smaller and larger fragments, respectively. To explain composite and single 



power law distributions, we can however use the models in |IU |. 

To conclude we introduced statistical models with a fixed probability for 
the breaking the largest particle in the system. All other particles of sizes 
smaller than the largest size, were divided into two different regions with 
different rates of fragmentation in each region. In the models therefore the 
largest fragments are broken with very large probability. This fact together 
with the form of the fragmentation rate as a function of the fragment size 
in each region, give a distribution of power law or log-normal form in that 
region. Our models therefore give a theoretical explanation for the transition 
of distribution from the scaling to log-normal behaviour, which had been seen 
before in the experiments. 
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